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^ 1 Abstract. The number density and clustering properties of Lyman-break galax- 
£3 . ies (LBGs) observed at redshift z ~ 3 are best explained by assuming that they 
are associated with the most massive haloes at z ~ 3 predicted in hierarchical 
models of structure formation. In this paper we study, under the same assump- 



tion, how star formation and chemical enrichment may have proceeded in the 
LBG population. A consistent model, in which the amount of cold gas available 
for star formation must be regulated, is suggested. It is found that gas cool- 
ing in dark haloes provides a natural regulation process. In this model, the star 
formation rate in an LBG host halo is roughly constant over about 1 Gyr. The 
predicted star formation rates and effective radii are consistent with observations. 
The metallicity of the gas associated with an LBG is roughly equal to the chem- 
ical yield, or about the order of 1Z@ for a Salpeter IMF. The contribution to the 
total metals of LBGs is roughly consistent with that obtained from the observed 
cosmic star formation history. The model predicts a marked radial metallicity 
gradient in a galaxy, with the gas in the outer region having much lower metallic- 
ity. As a result, the metallicities for the damped Lyman-alpha absorption systems 
expected from the LBG population are low. Since LBG halos are filled with hot 
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gas in this model, their contributions to the soft X-ray background and to the 
UV ionization background are calculated and discussed. 



Key words: galaxies: LBGs - galaxies: formation - galaxies: star formation - 
galaxies: chemical evolution 

1. Introduction 

The Lyman-break technique (e.g. Steidel, Pettini & Hamilton 1995) has now been 
proved very successful in finding large numbers of star forming galaxies at redshift 
z ~ 3 (e.g. Steidel et al. 1996, 1999b). The observed number density and clus- 
tering properties of Lyman-break galaxies (hereafter LBGs, Steidel et al. 1998; 
Giavalisco et al. 1998; Adelberger et al. 1998) are best explained by assuming 
that they are associated with the most massive haloes at z ~ 3 predicted in hi- 
erarchical models of structure formation (Mo & Fukugita 1996; Baugh, Cole & 
Frenk 1998; Mo, Mao & White 1998b; Coles, et al. 1998; Governato, et al. 1998; 
Jing 1998; Jing & Suto 1998; Katz, et al. 1998; Kauffmann, et al. 1998; Moscar- 
dini, et al. 1998; Peacock, et al. 1998; Wechsler, et al. 1998). This assumption 
provides a framework for predicting a variety of other observations for the LBG 
population. Steidel et al. (1999b and references therein) gave a good summary of 
recent studies on this population including the luminosity functions, luminosity 
densities, color distribution, star formation rates, clustering properties, and the 
differential evolution. 

Assuming that LBGs form when gas in dark haloes settles into rotationally 
supported discs or, in the case where the angular momentum of the gas is small, 
settles at the self-gravitating radius, Mo, Mao & White (1998b) predict sizes, 
kinematics and star formation rates and halo masses for LBGs, and find that the 
model predictions are consistent with the current (rather limited) observational 
data; Steidel et al. (1999a) suggest that the total integrated UV luminosity den- 
sities of LBGs are quite similar between redshift 3 and 4 although the slope of 
their luminosity function might have a large change in the faint-end. 

Furthermore, Steidel et al. (1999b) suggest that a "typical" LBG have a star 
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scale be the order of lGyr based on their values of E(B-V) as pointed out by Pet- 
tini et al. (1997b) after adopting the reddening law of Calzetti (1997). Recently, 
Friaca & Terlevich (1999) use their chemodynamical model to propose that an 
early stage (the first Gyr) of intense star formation in the evolution of massive 
spheroids could be identified as LBGs. 

However, Sawicki & Yee (1998) argued that LBGs could be very young stellar 
populations with the age less than 0.2Gyr based on the broadband optical and 
IR spectral energy distributions. This is also supported by the work of Ouchi & 
Yamada (1999) based on the expected sub-mm emission and dust properties. It 
is worthy of noting that the assumptions about the intrinsic LBG spectral shape 
and the reddening curve play important roles in these results. 

In this paper, we study how star formation and chemical enrichment may have 
proceeded in the LBG population. As we will demonstrate in Section 2, the 
observed star formation rate at z ~ 3 requires a self-regulating process to keep 
the gas supply for a sufficiently long time. We will show (in Section 2) that such 
a process can be achieved by the balance between the energy feedback from star 
formation and gas cooling. Model predictions for the LBG population and further 
discussions about the results are presented in Section 3, a brief summary is given 
in Section 4. 

As an illustration, we show theoretical results for a CDM model with cosmo- 
logical density parameter VLq = 0.3, cosmological constant = 0.7. The power 
spectrum is assumed to be that given in Bardeen et al. (1986), with shape pa- 
rameter T = 0.2 and with normalization a$ = 1.0. We denote the mass fraction 
in baryons by /b = ^b/^o, where f^B is the cosmic baryonic density parameter. 
According to the cosmic nucleosynthesis, the currently favoured value of f^B is 
S7b ~ 0.019/i~ 2 (Buries & Tytler 1998), where h is the present Hubble constant 
in units of 100 kms _1 Mpc _1 , and so /b ~ 0.063/i~ 2 . Whenever a numerical value 
of h is needed, we take h = 0.7. At the same time, we define parameter i* as the 
time scale for star formation in the LBG population throughout the paper. 



2. Models 
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2.1. Galaxy Formation 

In this paper, we use the galaxy formation scenario described in Mo, Mao & White 
(1998a, hereafter MMWa) to model the LBG population. In this scenario, central 
galaxies are assumed to form in dark matter haloes when collapse of protogalactic 
gas is halted either by its angular momentum, or by fragmentation as it becomes 
self-gravitating (see Mo, Mao & White 1998b, hereafter MMWb, for details). As 
described in MMWb, the observed properties of LBGs can be well reproduced if 
they are assumed to be the central galaxies formed in the most massive haloes 
with relatively small spins at z ~ 3. As in MMWb, we assume that gas in a dark 
halo initially settles into a disk with exponential surface density profile. 

When the collapsing gas is arrested by its spin, the central gas surface density 
and the scale length of an exponential disk are 
2 ( m d \ ( A \ ~ 2 / V c 



S » 38O/iM pc 
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where ma is the fraction of halo mass that settles into the disk, V c is the circular 
velocity of the halo, A is the dimensionless spin parameter, H(z) is the Hubble 
constant at redshift z and Hq is its present value (see MMWa for details). Since 
H(z) increases with z, for a given V c disks are less massive and smaller but have 
a higher surface density at higher redshift. When A is low and m d is high, the 
collapsing gas will become self-gravitating and fragment to form stars before it 
settles into a rotationally supported disk. In this case, we will take an effective 
spin A oc nid m calculating So and R^. 

We take the empirical law (Kennicutt 1998) of star formation rate (SFR) to 
model the star formation in high-redshift disks which is 

Es ™ = (i^)' Ms! "'" pc "' (3) 

where 

a = 2.5xl0- 10 , b = lA (4) 
respectively. Here SgpR is the SFR per unit area and S gas is the gas surface 
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formation rate and cold gas density over large areas on spiral disks and over 
starburst regions (Kennicutt 1998). We will apply this law differentially on a 
disk and also take into account the Toomre instability criterion of star formation 
(Toomre 1964; see also Binney & Tremaine 1987). 

For a given cosmogonic model, the mass function for dark matter haloes at red- 
shift z can be estimated from the Press-Schechter formalism (Press & Schechter 



where 5 c (z) = <5 C (0)(1 + z)g(0)/g(z) with g(z) being the linear growth factor at z 
and 5 C (0) ss 1.686, A(i?) is the linear rms mass fluctuation in top-hat windows 
of radius R which is related to the halo mass M by M = (47r/3)p i? 3 , with p 
being the mean mass density of the universe at z = 0. The halo mass M is related 
to halo circular velocity V c by M = /[10GH(z)]. A detailed description of the 
PS formalism and the related cosmogonic issues can be found in the Appendix 
of MMWa. 

From the Press-Schechter formalism and the A-distribution which is a log- 
normal function with mean In A = In 0.05 and dispersion <7i n A = 0.5 (see equation 
[15] in MMWa), we can generate Monte Carlo samples of the halo distributions 
in the V c -X plane at a given redshift and, using the star formation law outlined 
above, assign a star formation rate to each halo. As in MMWb, we select LBGs as 
the galaxies with the highest star formation rate, so that the comoving number 
density for LBGs is equal to the observed value, A^lbg = 2.4 x 10~ 3 /i 3 Mpc~ 3 
for the assumed cosmology at z = 3, as given in Adelberger et al. (1998). Here 
it is worth noting that the model selection of LBGs we adopted is without the 
dust extinction being considered. This implies that the contribution of the dust 
is assumed to be uniform. But in fact, it could be very different from galaxies to 
galaxies. So, our selection of LBGs may not have one-to-one correspondence with 
the observed LBGs (Baugh et al. 1999), but the selection should be correct on 
average. 



1974): 
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2.2. Cooling-Regulated Star Formation 

What regulates the amount of star-forming gas in a dark halo? In the standard 
hierarchical scenario of galaxy formation (e.g. White & Rees 1978; White & Frenk 
1991, hereafter WF), gas in a dark matter halo is assumed to be shock heated to 
the virial temperature, 

T = 2.24 x 10 6 K ( — -] , (6) 

V 250km s-V V ; 

as the halo collapses and virializes. The hot gas then cools and settles into the 

halo centre to form stars. As suggested in WF, the amount of cold gas available 

for star formation in a dark halo is either limited by gas infall or by gas cooling, 

depending on the mass of the halo. For the massive haloes (V c ?Z 200km s -1 ) we 

are interested here, gas cooling rate is smaller than gas-infall rate, and the supply 

of star- forming gas is limited by gas cooling (see WF for details). It is therefore 

likely that gas cooling is the main process that constantly regulates the SFR in 

LBGs. 

To have a quantitative assessment, let us compare different rates involved in 
the problem. Using equations (l)-(4) we can write the SFR as 

M *-—& 2.33x10/. [ 25Qkms - 1 ) [~^-_ 

where ma is the current gas content of the disk. The rate at which gas is consumed 
by star formation is therefore 

M SFR = (1 - R T )M„ (8) 

where R r is the returned fraction of stellar mass into the ISM; we take R T = 0.3 
for a Salpeter IMF (e.g. Madau et al. 1998). According to WF, the heating rate 
due to supernova explosions under the approximation of instantaneous recycling 
can be written as 
dE 

— = e M*(700kms- 1 ) 2 , (9) 

where eo is an efficiency parameter which is still very uncertain. We take it to be 
0.02 as in WF. The rate at which gas is heated up (to the virial temperature) is 
therefore 
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which is the same form as equation (9) of Kauffmann (1996; see also Somerville 
1997). At z = 3 and for the cosmology considered here, this rate can be written 

as 



M heat » 29.2/i 



-0.6 f m d \ 14 / A \-° V ^ \ U r^ 1_M 



0.05 7 V 0.05 7 V 250km s 



An 



Moyr- 1 ^!!) 



Comparing this equation with equations (7) and (8), we can find that the rate 
for gas consumption due to star formation is much larger than the rate of gas 
heating for LBG halos. Because LBGs are hosted by massive halos which have 
large circular velocities V c , the halos are cooling dominated which is confirmed 
during the detailed calculation below. Following WF we define a mass cooling 
rate by 

M coo i = 47rp gas (r coo i)r^ ool ^p, (12) 

where r coo i is the cooling radius and p gas is the density profile of the hot gas in the 
halo. For simplicity, we assume that p gas {r) = /bK? /(4irGr 2 ), and we define r coo i 
to be the radius at which the cooling time is equal to the age of the universe, 
which is similar to the time interval between major mergers of haloes (Lacey 
& Cole 1994). The density distribution of the halo mass here is assumed to be 
isothermal. However, it is the NFW profile (Navarro, Frenk & White, 1997) in 
MMWb. Because the difference of the resulted cooling rates between these two 
different choices of density profiles is small (Zhao et al, 1999), and the major goal 
here is to show whether or not the cooling-regulated star formation can be valid, 
the adoption of isothermal profile will not influence the final result very much. 

Under this definition, gas within the cooling radius can cool effectively before 
the halo merges into a larger system where it may be heated up to the new 
virial temperature if it is not converted into stars. Using the cooling function 
given by Binney & Tremaine (1987) where cooling function A « 10 _23 ergs _1 cm 3 
in the range of 5 x 10 5 K ^ T ^ 2 x 10 7 K (and assuming gas with primordial 
composition), the mass cooling rate can then be written as 

V r Wf R \ 3/2 



M, 



cool 



^ (250^) (M) M ^ _1 - <13) 



If M* is smaller than M coo \, then cold gas will accumulate in the halo centre 
and lead to higher star formation rate. If, on the other hand, M± > M coo \, the 
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leading to a lower star formation rate. We therefore assume that there is a rough 
balance among these three rates: 

M cool » M hcat + (1 - R T ) M*. (14) 

It should be noted that the cooling-regulated star formation process is only a 
reasonable hypothesis, and the real situation must be much more complicated. 
For example, during a major merger of galactic haloes, the amount of gas that 
can cool must be much larger than that given by the cooling argument, and the 
star formation may be in a short burst (e.g. Mihos & Hernquist 1996). However, 
such bursts are not expected to dominate the observed LBG population, because 
of their brief lifetimes. Thus, star formation rates in the majority of LBGs are 
expected to be regulated by equation (14) on average. As shown in MMWb, to 
match the observed number density of LBGs, the median value of V c is about 
300 km s -1 in the present cosmogony. The typical star formation rate is of the 
order lOOM yr _1 . This is not very different from the observed star formation 
rates, albeit dust distinction in the observations may be difficult to quantify. 

Figure 1 shows the value of ma required by the balance condition equation (14) 
as a function of halo circular velocity, assuming that /b = 0.1 and the left hand 
side exactly equals to the right hand ones in equation (14). Results are shown for 
two choices of spin parameters, A = 0.035 and 0.08, corresponding to the 50 and 
90 percent points of the A distribution for the LBG population (MMWb). As one 
can see, for the majority of LBG hosts, gas cooling indeed regulates the values 
of ?rtd to the range from 0.02 to 0.04. So, we can reasonably choose = 0.03 
for the LBG population as MMWb did. Since the cooling time is approximately 
the age of the universe at z ~ 3, cooling regulation ensures that star formation 
at the predicted rate can last over a large portion of a Hubble time. 

3. MODEL PREDICTIONS FOR THE LBG POPULATION 

Since the cooling regulation discussed above gives specific predictions of how star 
formation may have proceeded in LBGs, here we use this model to predict the 
properties of the LBG population. The condition in equation (14) implies that 
the star formation rate in a disk is equal to the rate of gas infall (due to a balance 
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Fig. 1. The value of required by the balance condition equation (14) as a 
function of halo circular velocity V c at z = 3 for A = 0.035 and A = 0.08, 
assuming /b = 0.1 (see text). 

host halo is described by the standard chemical evolution model with infall rate 
equal to star formation rate, i.e., the new infalling gas to the disk distributed 
radially in an exponential form with the scale length of Rd/b ~ 0.7-Rd, and the 
reheated gas removed decreases with the increasing radius due to the decreasing 
SFR. Under the instantaneous recycling approximation (Tinsley 1980), the gas 
metallicity Z is given by 

Z = y(l-e- v ) + Zi, v =^-l, (15) 

where Z\ is the initial metallicity of the infalling gas, y is the stellar chemical 
yield, S gas is the gas surface density (which is kept constant by gas infall) and 
E to t is the total mass surface density, which increases as star formation proceeds: 

= (1 - ^)S SF R. (16) 

Here the enrichment of the halo hot gases is not taken into account because the 
amount of metals heated up to the halos by SNs is relatively smaller than that 
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Fig. 2. Predicted SFR as a function of V c and A in the cooling-regulated model, 
(a) SFR vs V c , for A =0.03, 0.05 and 0.1 (from top to bottom), (b) SFR vs A, for 
V c =300, 200 and lOOkm/s (from top to bottom). 

3.1. Individual Objects 

Figure 2 shows the star formation rate as a function of halo circular velocity 
V c and spin parameter A. As expected, the predicted SFR increases with V c but 
decreases with A . As we can see from the figure, if we define systems with 
SFR ^ 4OM yr _1 (which matches the SFRs for the observed LBG population) 
to be LBGs, the majority of their host haloes must have V c ^ 200 km s -1 which 
are cooling dominated. This result is the same as that obtained by MMWb based 
on the observed number density and clustering of LBGs. Thus, the star formation 
rate based on cooling argument is also consistent with the observed number den- 
sity and clustering. Because SFR is higher in a system with smaller A, the LBG 
population are biased towards haloes with small spins, but given its relatively 
narrow distribution, this bias is not very strong. 

The predicted metallicity gradients on individual disks are shown in Figure 3 
for two different choices of star formation time scale t* of 0.5Gyr and lGyr respec- 
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Fig. 3. The metallicity gradients for LBGs for different star formation time t* 
assuming that y = Z Q and Z\ = (see text). Full and dash lines show results for 
V c = SOOkms^ 1 and 150kms _1 , respectively. From top to bottom, A = 0.03 and 
0.1; (a) U =0.5Gyr; (b) U =lGyr 



easily compare with observations. The metallicity gradients are negative in all 
cases. When radius is measured in disk scale length, the predicted metallicity 
depends weakly on V c but strongly on A, and is higher for a longer star formation 
time. As one can see from equation (15), the largest metallicity in the model is 
Z = Zi + y. This metallicity can be achieved in the inner part of compact disks 
(with small A) when star formation time t+ ^> 1 Gyr. The metallicity drops by a 
factor of ~ 2 from its central value at R ~ 3i?d- 
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Fig. 4. The predicted metallicity distributions for LBG populations assuming 
that y = Z & and Z\ = in order to make the predictions easily compare with 
observations (see text). Results are shown for two star formation timescales £* = 
0.5 Gyr (dash) and t+ = 1 Gyr (solid), respectively (cf. equation (15)). 

3.2. LBG Population 

Since the distribution of haloes with respect to V c and A are known, we can 
generate Monte-Carlo samples of the halo distribution in the V c -\ plane at any 
given redshift. We can then use the galaxy formation model (MMWb) discussed 
above to transform the halo population into an LBG population based on LBGs 
with highest SFRs which is the same as that outlined in Sec. 2. 

We define the typical metallicity of a galaxy as the one at its effective radius. 
Figure 4 shows the distribution of this metallicity for two choices of the star 

r„ ... x: -t n tr /~i — 1 m — T..„x - „ „ _ _ . tt>: o 
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Fig. 5. The predicted effective-radius distribution for LBGs in the cooling- 
regulated scenario (solid), compared to the observed distribution (dash). 

in last section, we have assumed that y = Z & and Z\ = in order to make the 
predictions easily compare with observations. The median values of (Z — Zi)/y 
are 0.60 and 0.84 for t± = 0.5 Gyr and 1 Gyr, respectively. The sharp truncation 
at (Z — Zi)/y = 1 is due to the fact that this quantity has a maximum value of 
1 in the present chemical evolution model. It can be inferred form Figure 3 that 
the range in (Z — Zi)/y decreases with increasing star formation time. Thus, if 
gas infall lasts for a long enough time, the distribution in (Z — Z,j)/y will be very 
narrow near 1 and all LBGs will have metallicity Z = Zi + y. According to the 
works of Tinsley (1980) and Maeder (1992), the stellar yield y is the order of Zq 
for the Salpeter IMF. If we adopt a stellar yield y ~ 0.5Z Q and Zi = 0.01Z Q , 
and if LBGs are not short bursts (e.g. ^ 0.5 Gyr) then their metallicity will 
be Z ^ 0.2Z G which is similar to that proposed by Pettini (1999). 

The predicted distribution of effective radii for the LBG population is shown 
in Figure 5. The distribution is similar to that of MMWb. The predicted range is 
1.0 5s R c g 5s 5.0 fo _1 kpc with a median value of 2.5 /i _1 kpc. Note that the effective 
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Fig. 6. The predicted SFR distribution for LBGs in the cooling-regulated sce- 
nario. 

t+ and m&. The model prediction is in agreement with the observational results 
given by Pettini, et al. (1998), Lowenthal, et al. (1997) and Giavalisco et al. (1996) 
which are mentioned above. 

The predicted SFR distribution of LBGs also resembles the prediction of 
MMWb except for a slight difference with MMWb, which is shown in Figure 
6. The median values are 18OM yr _1 for the model and spans from 100 to 
5OOM yr _1 . To compare with observations, we have to take into account the 
effect of dust. If we apply an average factor of 3 in dust extinction, then the pre- 
dictions closely match the values derived from infrared observations by Pettini, 
et al. (1998) although there might exist rare LBGs with very high SFR. 

3.3. Contribution To The Soft X-ray and UV Background 

Since the virial temperature of LBG haloes are quite high, in the range of 10 6 — 
10 7 K, significant soft X-ray and hard UV photons may be emitted as the halo 
hot gas cools. It is therefore interesting to examine whether the LBG population 
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The dominant cooling mechanism for hot gas with temperature ^ 10 6 K is the 
thermal bremsstrahhmg. The bremsstrahlung emissivity is given by (e.g., Peebles 
1993) 

j„ = 5.4 x 10- 39 n2T- 1 / 2 e -W^ ergcm -3 s -i gter -i Hz -i > ( 17 ) 

where n e (in cm" 3 ) is the electron density and T (in K) is the temperature given 
by equation (6). The total power emitted per unit volume is 

J = 1.42 x 1(T 27 T 1 / 2 n 2 erg cm" 3 s _1 . (18) 

We write the total luminosity Lb in thermal bremsstrahlung as 

L h = (5M coo xVl (19) 

and we take (3 = 2.5 here as WF so that L\> is equal to the initial thermal energy in 
the cooling gas. Note that the value of (5 is quite uncertain because it depends on 
the detail density and temperature profiles of the hot gas. Substituting equation 
(13) into the above equation, we obtain the total soft X-ray luminosity for an 
LBG 

yMj: ( V c V (fa \ 3/2 _„ „-l 



f -J_r 

/soft " kTJo.. 



r ~ /kT dE (21) 



L„ W )«4.1xlO«,^^j ^) ergs-, (20) 
where 

r2(i+«) „ 

5(1 + 2) 

is the fraction of total energy that falls into the ROSAT soft X-ray (0.5-2 keV) 
band. The contribution of the LBG population to the soft X-ray background is 
then 

Psx = / / dWcom^f « 5.7 x 10- 8 (^) 3/2 erg s^W 2 , (22) 

where n(z) is the comoving number density of LBG haloes as a function of redshift 
z, dV com is the differential comoving volume from z to z + dz and di, is the 
luminosity distance. The integrate for V c is to sum up all selected LBGs with 
V c based on their highest SFRs. We have integrated over redshift range from 3 
to 4 where the number density of LBGs is nearly a constant (Steidel 1998a,b). 
This contribution should be compared with the value derived from the ROSAT 
observations (Hasinger et al. 1998) in the 0.5-2 keV band 
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As we can see, the soft X-ray contribution from LBGs could be a substantial 
fraction (about 20%) of the total soft X-ray background. 

Similarly we can calculate the contribution of LBGs to the UV background at 
z = 3. We evaluate the UV background at 4 Ryd (lRyd=13.6 eV) using nearly 
identical procedures, we find that 



which is much smaller than the UV background from AGNs, i4R y d ~ 
10 _22 ergs _1 cm _2 Hz _1 ster~ 1 (e.g. Miralda-Escude & Ostriker 1990). 

3.4- Contribution to the Total Metals 

Based on the recent observational results of the cosmic star formation history, 
Pettini (1999) obtained a predicted total mass of metals produced at z = 2.5. 
After combining results of all contributors observed, he argued that there seems 
to exist a very serious "missing metal" problem, i.e., the predicted result is much 
higher than observed ones. So, it is interesting to evaluate the total metals pro- 
duced by LBGs in our model. 

According to the method we select LBGs to be the galaxies with highest SFR 
and our chemical evolution model mentioned in Sec. 3.2, we can calculate the 
total metal density produced by the LBG population at z = 3 based on their 
observed comoving number density which is TVlbg = 2.4 x 10~ 3 /i 3 Mpc~ 3 for 
the assumed cosmology (Adelberger et al. 1998). Defining that Qz is the metal 
density relative to the critical density, we get that fiz of LBGs are 0.19^6 x U 
and 0.2917b x V for star formation time of 0.5Gyr and lGyr respectively, where 
y is the stellar yield which is the same as above. Because the virial temperature 
of LBG halos are very high, a significant fraction of the metal should be in hot 
phase. Comparing our results with that estimated by Pettini (1999) which is 
0.08^6 x Z & (the cosmogony has been taken into account), we find that there is 
no "missing metal" problem in our model. 

3.5. LBGs and Damped Lyman- Alpha Systems 

Damped Lyman-alpha systems (DLSs) are another population of objects that 
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to their high neutral HI column density (> 10 20 3 cm -2 ), and are believed to 
be either high-redshift thick disk galaxies (Prochaska & Wolfe 1998) or merg- 
ing protogalactic clumps (Haehnelt, Steimetz & Rauch 1998). In either case, to 
match the observed abundance of DLSs, most DLSs should have circular velocity 
between 50kms -1 to 200 km s -1 , much smaller than the median circular veloc- 
ity of LBGs (~ 300 kms -1 ). Based on the PS formalism (equation (5)) and disk 
galaxy formation scenario suggested by MMWa (equations (1) and (2)), we can 
estimate with the random inclination being taken into account, that the fraction 
of absorbing cross-sections contributed by LBGs amounts to only about 5% of 
the total absorption cross-section assumed LBGs with highest SFRs. This means 
that only a very small fraction of DLSs can be identified as LBGs. 

The physical connection between LBGs and DLSs is still unclear, although 
the recent observation of Moller & Warren (1998) using HST indicates that 
some DLSs could be associated with LBGs. In Figure 7, we show the predicted 
metallicity distribution for the subset of DLSs which can be observed as LBGs. 
aGain, we have assumed that y = Z & and Z\ = to let the predictions more 
easily compare to observations. As can be seen, the DLSs generally have lower 
metallicity than LBGs, because they are biased towards the outer region of the 
host galaxies, where the star formation activity is reduced. Notice, however, that 
the metallicity of these DLSs could still be higher than most DLSs at the same 
redshift, which typically have metallicity of 0.1Z® (Pettini, et al 1997a). 

4. SUMMARY 

In this paper, we have examined the star formation and chemical enrichment 
in Lyman break galaxies, assuming them to be the central galaxies of massive 
haloes at z ~ 3 and using simple chemical evolution models. We found that gas 
cooling in dark haloes provides a natural process which regulates the amount 
of star forming gas. The predicted star formation rates and effective radii are 
consistent with observations. The metallicity of the gas associated with an LBG 
is roughly equal to the chemical yield, or about the order of 1Z® for a Salpeter 
IMF. Because of the relatively long star-formation time, the colours of these 
galaxies should be redder than that of short starbursts. It is not clear whether 



18 



Chenggang Shu: Star Formation and Chemical Evolution of LBGs 




Fig. 7. The predicted metallicity distribution for the DLSs expected from the 
LBG population (see text). Results are shown for two star formation timescales 
i* = 0.5 Gyr (dash) and = 1 Gyr (solid), respectively. 

the interpretation of the observational data depends strongly on the adopted dust 
reddening. Stringent constraint can be obtained when full spectral information 
of the LBG population is carefully analyzed. 

The model predicts a marked radial metallicity gradient in an LBG, with the 
gas in the outer region having lower metallicity. As a result, the metallicities for 
the damped Lyman-alpha absorption systems expected from the LBG popula- 
tion are lower than those for the LBGs themselves, although high metallicity is 
expected for a small number of sightlines going through the central regions of an 
LBG. At the same time, our modeled contribution to the total metal is roughly 
consistent with that obtained from the observed cosmic star formation history, 
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be more than half of the metals to be in the hot phase. Finally, a prediction of 
our model is that LBG haloes are filled with hot gas. As a result, these galax- 
ies may have a non-negligible contribution to the soft X-ray background. The 
contribution of LBGs to the ionizing UV background is found to be small. 

There are two basic assumptions in our work. One is that the LBG population 
is one-to-one associated with the most massive halos which are generated from 
the PS formalism, as done by MMWb; another is that the timescale of star 
formation for LBG population is assumed to be the order of lGyr, which is 
suggested by Steidel et al. (1999a,b, 1995). However, Baugh et al. (1999) recently 
argue that the prediction of the clustering properties of LBGs based on this 
first simple assumption will be discrepancy with the results of more detailed 
semi-analytic models. Still, the second will lead to difficulty in reproducing the 
redshift evolution of bright galaxies (Kolatt et al. 1999). More detailed modelling 
done by Somerville (1997) suggest the collisional starbursts could be expected to 
be an important effect in understanding the LBGs. So, further observations are 
required to investigate the intrinsic properties of LBGs. 
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